Quantum-Dot Cellular Automata using Buried Dopants 
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The use of buried dopants to construct quantum-dot cellular automata is investigated as an 
alternative to conventional electronic devices for information transport and elementary computation. 
This provides a limit in terms of miniaturisation for this type of system as each potential well is 
formed by a single dopant atom. As an example, phosphorous donors in silicon are found to have 
good energy level separation with incoherent switching times of the order of microseconds. However, 
we also illustrate the possibility of ultra-fast quantum coherent switching via adiabatic evolution. 
The switching speeds are numerically calculated and found to be 10's of picoseconds or less for a 
single cell. The effect of decoherence is also simulated in the form of a dephasing process and limits 
are estimated for operation with finite dephasing. The advantages and limitations of this scheme 
over the more conventional quantum-dot based scheme are discussed. The use of a buried donor 
cellular automata system is also discussed as an architecture for testing several aspects of buried 
donor based quantum computing schemes. 

PACS numbers: 85.35.Gv,85.35.p,03.67.Lx 



I. INTRODUCTION 

The use of quantum systems to perform computing 
tasks is an area of continued interest, both in the con- 
text of fully coherent quantum computers and the more 
general area of nano-computing. Recently there has been 
much interest in producing an experimental analogue of 
cellular automata at the micro or even nanometer scaled. 
This could provide an alternative architecture with which 
to build standard logic gates and information channels, 
which have low dissipation-, and fewer control gates com- 
pared to conventional transistor based logic. Significant 
progress has been made in this area with the experimen- 
tal demonstration of Quantum-Dot Cellular Automata 
(QDCA) cells constructed from Aluminium quantum- 
dots^ and GaAs/AlGaAs heterostructuresA. Recently, 
the operation of the functional components of a QDCA 
logic gate has been demonstrated experimentally 7 . This 
type of system has also been investigated as a possible ar- 
chitecture for quantum computing using QDCA qubits* 
and more recently as a candidate for a decoherence free 
subspaca^. 

In this paper we explore an alternative QDCA archi- 
tecture using buried dopants in semiconductors, as this 
provides a very strongly confined potential and well char- 
acterised energy levels. A recent proposal describes a 
charge-based qubit for quantum computing using phos- 
phorous donors in siliconifliii. We show that this treat- 
ment can be applied to a system of dopants arranged 
in the layout of a QDCA. We refer to these structures 
as Buried Dopant Cellular Automata or BDCA. Such 
a device could be fabricated by either direct atomic 
placementiSii^ or ion-implantationiiii4. Numerical es- 
timates are given for the case of phosphorous donors 
in silicon but the concepts are generally applicable to 
other dopants. This is of particular interest given the 
recent advances in single dopant placement and cluster 



based charge transfer experiments using phosphorous in 
silicon^. 

Conventional QDCA rely on incoherent evolution (gov- 
erned by the T\ relaxation time) to mediate transitions 
between the logical states. For phosphorous donors in 
silicon we estimate this relaxation time and find it to 
be of the order of microseconds. We also investigate 
an alternative ultra-fast (picosecond) switching mecha- 
nism, namely coherent evolution between defined sys- 
tem eigenstates. This approach is central to the use 
of buried donors and is also applicable to coherent tun- 
nelling between quantum dots or superconducting sys- 
tems. This constitutes an alternative evolution mecha- 
nism for QDCA schemes where coherence can be main- 
tained long enough for the cell to be switched from one 
classical state to another without the need for long co- 
herence times which are typically required for quantum 
computing applications. The relevant time-scales and ap- 
propriate pulsing sequences with and without dephasing 
are discussed and the scaling behaviour of the system is 
investigated. 



II. QUANTUM-DOT CELLULAR AUTOMATA 

The simplest QDCA is a cell composed of four 
quantum-dots containing two mobile electrons which can 
move between the dots via tunnel junctions. The elec- 
trons tend to occupy diagonally opposite sites to min- 
imise the energy due to the Coulombic interaction. These 
two ground (or computational) states are labelled zero 
and one, see Fig. ^ where 'e' indicates the position 
of the electrons. The next highest energetic states are 
non-computational states and ideally are only transiently 
populated during correct operation. 

If two dots are placed next to each other, one cell in- 
fluences the state of the other cell via capacitative cou- 
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FIG. 1: Two possible states for a basic QDCA cell where the 
and 1 states constitute the ground or 'computational' states 
and 'e' labels the position of the electrons. 
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FIG. 2: Layout for a QDCA wire (a) and inverter (b) which 
demonstrate information transfer and binary inversion respec- 
tively after Tougaw and Len*-. 



pling. In an array of cells, when the first cell is switched 
from one computational state to the other, the rest of the 
chain relaxes to minimise the energy of the total system. 
The result of this relaxation is to transfer the state in- 
formation of the initial cell along the chain without net 
electron flow and minimal energy dissipation. The speed 
at which this switching occurs is governed by the incoher- 
ent tunnelling rate of the junctions, the inverse of which 
is referred to as the Tj or relaxation time. Classical in- 
formation processing can be performed in this scheme, 
as shown in Fig. [21 for a QDCA wire and inverter 1 . It is 
also possible to realise non-trivial classical computation, 
such as a full-adder, using this scheme^. 



III. BURIED DONORS AND THE 
HYDROGENIC APPROXIMATION 

The use of buried donors in a semiconductor ma- 
trix has been discussed for charge-based quantum 
computingiSii^ii^. While the advantages of semiconduc- 
tor fabrication and gate control are well known, the fast 
dephasing and relaxation effects mean that charge-based 
quantum computing using buried donors is still techni- 
cally difficult. On the other hand the QDCA architec- 
ture is not as seriously affected by dephasing or relax- 
ation, as the system is always in the ground state when 
it is measured. We will investigate the possibility of us- 
ing this concept for building QDCA based systems using 
phosphorous donors in silicon (Si:P), which we term a 
BDCA. The basic idea is to construct an array of four 
ionised donors which contains two 'free' electrons, there- 



fore mimicking the layout of a conventional QDCA. This 
also represents a limit in terms of miniaturisation for this 
form of nano-computing as each potential well is created 
by only one donor atom. 

In order to provide a model for the Si:P donor system, 
we will use the effective mass or 'hydrogenic' approxima- 
tion in which the outer shell electron of a phosphorous 
donor in silicon can be treated as a hydrogen orbital with 
the energies and distances scaled appropriately. A more 
complete treatment of effective mass theory for shallow 
donors is given by Kohnit Eq. (Q) and @ give the scal- 
ing factors for the effective Bohr radius (a|j) and effective 
energy (E*) of the donor electron in terms of the effective 
mass of the donor electron and the dielectric constant of 
the substrate, 
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* m* 1 
E* = -E. 



(1) 



(2) 



The advantage of this approach is that solving a hydro- 
genic system with a small number of electrons is more 
tractable than a full electron calculation of the phospho- 
rous donor within a silicon lattice. These results also 
generalise to other shallow donor systems. We will con- 
centrate on phosphorous donors in silicon and therefore 
quote the appropriate energy levels, where the conver- 
sion is 13.6eV in a hydrogenic system is approximately 
equal to 20meV for Si:P. The effective Bohr radius of the 
phosphorous electron is approximately 3nm for an effec- 
tive electron mass m* = 0.2m e and a dielectric constant 
e = 11.7. 



IV. EFFECTIVE HAMILTONIAN 

To provide a convenient formalism, we construct an 
effective Hamiltonian using the pseudo-spin approach to 
describe the BDCA system. By defining each pair of 
phosphorous donors and their shared electron as a sin- 
gle pseudo-spin object we can define two states, top (T) 
and bottom (B), which specify the position of the elec- 
tron. Each BDCA cell then consists of a pair of these 
objects where the computational states are \TB) = |0) 
and \BT) = |1) respectively, as shown in Fig.[3{a). 

In this way, the labelling for the non-computational 
states is \TT) and \BB), as shown in Fig.[3Jb). Initially 
we will assume that the electrons cannot move laterally 
so each electron is bound to its particular donor pair, 
as indicated by the ellipses in Fig. [31 This would corre- 
spond to a situation where the direction of tunnelling is 
controlled by confining potentials or the geometry of the 
cell. A more complete justification for this assumption is 
given later. 

Once the position of the electrons is encoded using 
this pseudo-spin approach, an effective Hamiltonian can 
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FIG. 3: (a) The ground states of the buried donor BDCA cell 
where the positions of the electrons 'e' are designated by top 
(T) and bottom (B). These computational states are referred 
to as \TB) and \BT), and are assigned the logical values of 
and 1 respectively, (b) The excited or 'non-computational' 
states are labelled \TT) and \BB) respectively and correspond 
to the first excited state of the system. 



be developed using the Pauli spin representation, 
H cS =ei(Ti + e 2 (T2 + Aiof + A 2 fi2 

'2 I 

where the usual Pauli matrices arc defined as 



JxzCTi CT 2 + Jz^\^2 



(3) 
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The notation cr^ is used to refer to the matrix applied 
to the ith donor pair and the coefficients e, A and J are 
determined as follows. 

The true hydrogenic Hamiltoniani& of the two elec- 
tron/four donor system is given (in scaled atomic units) 
by 



H = — 



Vf - Vi - 2( 1 1 1 
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7*12 



(4) 



where is the separation between the ith electron and 
the jth donor (j = a, 6, c, d) and r\2 is the separation be- 
tween the electrons. Numerically evaluating this Hamil- 
tonian within the basis of four states (\TB), \BT), \TT) 
and \BB)) enables the elements of the matrix 



Hi. 



(5) 



to be found. The elements of Eq. JSJ are then equated to 
the coefficients in Eq. © to determine an effective Hamil- 
tonian. The basis of states is represented using a lin- 
ear combination of atomic orbitals (LCAO). We use the 
anti-symmetric spatial wavefunction for the H 2 molecule 
as our basis wavefunction. The spin-orbit coupling for 
donor electrons in silicon is known to be very small 19 




i — 1 — i — 1 — i — ' — i — ■ — i — ' — i — ■ — i — ' — r 



Numerical 
Calculation 
Electrostatic 
Approximation 



10 20 30 40 50 60 70 80 
Side Length R [nm] 



FIG. 4: Energy difference between the ground (E ss ) and first 
excited (-E ex ) states of a square BDCA cell made up of phos- 
phorous donors in silicon. The energy difference is computed 
for various donor separations (R) by numerically integrating 
the schrodinger equation. The points are full quantum me- 
chanical calculations using the LCAO approach and the solid 
line is the energy difference determined analytically from sim- 
ple electrostatic arguments, Eq. ©. 



and therefore spin effects may be neglected as the spin 
and charge degrees of freedom are assumed to be sepa- 
rable at all times. The elements of Eq. are then used 
to obtain estimates for the numerical coefficients in the 
effective hamiltonian (H c g) for a square BDCA cell of 
side length R. 

We estimate the energy difference (E cx — E gs ) between 
the excited (non-computational) and ground (computa- 
tional) states of one cell as a function of system size. 
This energy gap gives an estimate of the temperature 
at which the system must be operated to ensure that the 
non-computational states are not thermally excited. The 
energy difference for a range of separations is plotted in 
Fig. 0] The points are found using the LCAO approach 
while the line is an approximation found using electro- 
static arguments based on the geometry of the system 
and has the form 



(2 - V2)E*a\ 
R 



(6) 



where R is the side length of the BDCA cell in nm and 
a* B = 3nm, E* — 20meV are the effective Bohr radius 
and effective ground state energy respectively. This ap- 
proximation is found to be valid in the region where 
the electron wavefunctions do not strongly overlap with 
each other but deviates from the numerical results for 
R < lOnm. 

Assuming a Boltzmann distribution for the energy 
states at finite temperature, the occupation probability 
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(P ex ) of the excited states (\TT) and \BB)) is given by 

P ex «e~^fr, (7) 

where E cx — E gs = AE > is the energy difference be- 
tween the ground and excited states. For a separation of 
15nm, we require an operating temperature of < 3K to 
achieve > 99% occupation of the computational states. 
For an operating temperature of lOOmK, the occupancy 
of the computational states is approximately 100%. 

For the rest of the discussion, we will use the numeri- 
cal coefficients calculated for a square cell of side length 
15nm. As this is equivalent to a cell size of approximately 
5 Bohr radii, the electrons can be said to be well localised 
and overlap effects are not significant. In order to sim- 
plify the situation, we will also set the <r x , a x a x and 
a v a v type terms to zero and then reintroduce them later 
in a systematic fashion. This corresponds to a situation 
where a surface gate potential is used to control over- 
lap of the electron wavefunctions and therefore control 
the tunnelling rate between donors. For large cells sizes, 
this is a good approximation to the physical situation as 
the tunnelling rate without an applied barrier bias will 
be very low. For small cell sizes there will be wavefunc- 
tion overlap even without a barrier gate. In this case a 
confining potential can be used to localise the electrons 
and provide more control over the tunnelling character- 
istics. Given these approximations, many of the terms in 
Eq. @ are approximately zero and the only significant 
term for R = 15nm is J zz = 1.21 meV which is due to 
the Coulombic repulsion between the electrons. 



V. BDCA SWITCHING 

To study how a BDCA cell would switch from one 
computational state to another, we consider the effect 
of control gates and nearby cells. We propose a struc- 
ture where pairs of donors are positioned in a line with 
a surface 'barrier' gate constructed between them which 
is used to control the tunnelling rate between the pairs 
of donors. At the end of the chain, 'symmetry' gates al- 
low the system to be switched from one state to another, 
depending on the bias applied. Sensitive electrometers, 
such as single-electron transistors (SET)2£, are used to 
measure the position of the electron at the other end of 
the chain. A diagram of this concept is shown in Fig. [5] 

Most of the following section refers to a single BDCA 
cell made from two pairs of donors for simplicity, but the 
discussion applies equally to a long chain of donor pairs. 
The effect of the barrier (B) and symmetry gates (S) for 
this cell can be included by adding terms to the effective 
Hamiltonian, Eq. J2J, of the form 

E s (t)at-E B (t)(af + a%), (8) 

where the magnitude of the coefficients (E$ and Eb) are 
controlled by the voltages applied to the gates. 
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FIG. 5: Simplified layout of a BDCA chain, where the cell size 
is labelled R and the cell spacing (S) is the distance between 
the centre of neighbouring cells. The symmetry and barrier 
surface gates are used to perform switching and control the 
tunnelling rate respectively. The circles represent the position 
of the donors. The position of the electrons 'e' is shown for the 
ground state configuration which corresponds to some non- 
zero bias on the control gates with the labelled polarity. The 
position of the electrons are measured with single-electron 
transistors (SET). 



The symmetry gates localise the system in one of the 
two computational states, based on the gate's polarity. 
The voltages applied to each symmetry gate have equal 
magnitude but opposite sign to ensure a symmetrical ef- 
fect on the chain. 

The barrier gate controls the tunnelling rate between 
pairs of donors by repelling or attracting the electron 
clouds. If there is a large enough separation between 
donors, the barrier gate is required to allow tunnelling 
by reducing the potential barrier that the electron feels. 
This justifies our initial assumption that when the sep- 
aration is large enough tunnelling only occurs between 
donor pairs and not along the BDCA chain. The use 
of compensating gates could also be used to confine the 
electrons and therefore prevent tunnelling in unwanted 
directions. 

The effect of these surface gates has been modelled as 
pure a x and <r 2 terms in the Hamiltonian due to the sym- 
metry of the system. We estimate the surface gate volt- 
ages to be 100 — lOOOmV depending on the donor depth 
and the presence of an oxide barrier layer, based on es- 
timated gate voltages for charge-qubitsi2i2i. Assuming 
arbitrary high precision in donor placement, the barrier 
gate has an equal effect on donors either side of the bar- 
rier and can be considered a pure a x gate. The barrier 
gate would also induce o x o x and a y a v style coupling but 
this is expected to be small compared to the Jzz and pure 
a x coupling. In addition, the barrier gate can also have 
a negative bias applied to improve the localisation of the 
computational states during readout, though this is not 
directly modelled here. While the symmetry gates would 
not be pure <j z (having some residual a x effect due to 
wave-function overlap) this could also be compensated 
for by using the barrier gate or additional compensation 
gates. 

In Fig. the eigenspectrum for a single BDCA cell is 
plotted as a function of the symmetry gate potential (Eg) 
for zero barrier gate potential. The computational states 
(\TB) and \BT)) are localised even for a very small sym- 
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FIG. 6: Eigenspectrum for a four donor cell as the symmetry 
potential (Es) is swept from -2 to 2 (meV) with no applied 
barrier potential (Eb = 0). 



metry potential. Also note that the two computational 
states are degenerate when there is no potential differ- 
ence applied to the control gates, as expected. 



VI. INCOHERENT SWITCHING 

If we allow the system to evolve via incoherent relax- 
ation (in direct analogy with the quantum-dot schemes), 
the transition from a high to a low energy eigenstate of 
the system is mediated by phonons in the lattice. For 
the moment wc will ignore the effect of the barrier gate 
and assume the electrons are well localised. This situa- 
tion is valid over time scales greater than the dephasing 
time of the donor wavefunction or where the barrier gate 
has a negative bias producing a high potential barrier. 
To obtain an estimate for this rate, we assume that the 
electrons within a BDCA cell can relax independently 
and we therefore use a similar approach to that used to 
estimate the relaxation rate for charge-qubits based on 
buried donors 2 ^ 3 .. Following the approach used by Bar- 
rett and MilburnSS, and Bockelmann and Bastard^i, we 
write the relaxation rate due to thermal phonons, 



r p h — 



64D 2 qf f [n B (E,T ph ) + a][l - smc(q if R)] 
nphc 2 [(q lf a B ) 2 + A^ 



(9) 



where a = 1 for emission and for absorption of a 
phonon, qi$ is the wavenumber of a phonon with a mag- 
nitude equal to the energy difference between the states 
(I Si/ 1 = E/hc s ), R is the separation between the donors 
and ub(E, T p h) = [exp (£/fcflTph) — is the Bose oc- 
cupation function for a bath of phonons at temperature 
Tph- We have ignored effects due to coherent tunnelling 
and used the following parameters^ for Si:P where D = 
3.3eV is the deformation potential, p = 2329kgm~ 3 is 
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FIG. 7: Energy level diagram for a single BDCA cell in the 
presence of a well localised neighbouring cell, where the cell 
size R = 15nm and there is a cell spacing S = 30nm. The 
direct transition (\BT) — > \TB)) is suppressed as the inter- 
action is phonon mediated and must therefore proceed via 
single-electron transitions. The two (first order) decay paths 
from one computational state to the other are illustrated with 
their associated transition rates for absorption (A) and emis- 
sion (E). 



the density of silicon, c s = 9.0 x 10 3 ms _1 is the speed 
of sound in silicon and a B = 3nm is the effective Bohr 
radius of the donor electron. 

To estimate the incoherent switching time, we calcu- 
late the energy levels for a BDCA cell (the target cell) 
assuming that a neighbouring cell is well localised. In the 
incoherent limit, we can assume that each electron is well 
localised and so the effect of the neighbour cell is to lift 
the degeneracy of the computational states of the target 
cell. For a cell of side length R = 15nm and cell spacing 
S = 30nm, a well localised neighbour induces an energy 
splitting of 1.64meV between the computational states 
of the target cell. This is calculated using the difference 
in electrostatic repulsion between the target cell and its 
neighbour. The splitting caused by the neighbouring cell 
can be modelled as a bias on the symmetry gate of the 
target cell. In this case the equivalent symmetry gate 
bias is Es = 0.82meV. The resulting energy levels can 
be read from Fig. [fj] Using these energy levels we can 
calculate the relaxation rate from the first excited state 
to the ground state (\BT) — > \TB)) and estimate the 
switching time of the system. Relaxation in this system 
is phonon mediated and therefore acts on each electron 
separately. This means the direct transition from state 
\BT) to state \TB) is suppressed and instead the sys- 
tem must relax via a co-tunnelling (two electron) process 
which requires absorption and emission of phonons to 
reach the ground state. The two possible (first order) 
decay paths are illustrated in Fig. [7| for the energy levels 
corresponding to a symmetry potential of 0.82meV. 

To estimate the total relaxation rate, we add the co- 
tunnelling ratesSS for the two possible paths, 



Total 



r (i)p(i) s 

A E 



r (2) r (2W 
A E 



\E 



\TT) 



E 



\BT)\ 



E 



\BB) 



E 



\BT)\ 



(10) 



where T^'e are defined in Fig. and Bw is the energy 
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FIG. 8: Incoherent switching time calculated for a range of 
cell sizes (R) with the spacing between the neighbouring cells 
given by the cell centre-to-centre distance S = 2R. Higher 
operating temperatures result in faster switching times but 
also result in higher excited state populations, reducing the 
overall readout fidelity. 



of the |fc)th state. For these energy levels and an oper- 
ating temperature of 3K, the calculated relaxation rate 
TTotai = 1.1MHz, which gives a switching time of 0.9/is. 
This is almost two orders of magnitude slower than the 
estimated maximum switching rates of 90MHz for Al/Al- 
oxide QDCA structures^ at 70mK. This is to be ex- 
pected as there are no defined tunnel junctions in the 
buried donor case. The switching time is shown in Fig. El 
for a range of operating temperatures and cell sizes (R) 
with the spacing between neighbouring cells given by the 
cell centre-to-centre distance S = 2R. While the switch- 
ing rate does vary with cell size, the temperature effects 
dominate in this regime as the system requires enough 
thermal energy to mediate the two-electron transition. 
At higher temperatures, faster switching is expected but 
the occupation of the ground state is reduced, according 
to Eq. Q, resulting in an overall loss of fidelity. 

While the switching time is slow compared to modern 
microelectronics, it does demonstrate that if the system 
is initially setup in some excited state it will decay to the 
ground state. This process could be used to intialise the 
system in a known state. 



VII. COHERENT SWITCHING 

An alternative to the incoherent switching is to use 
the high tunnelling rates of charge-based quantum- 
computing schemes to perform coherent switching of the 
BDCA chain. As we have relatively strong coupling be- 
tween electrons in this system, we can consider adia- 
batic evolution as a mechanism to switch from one com- 
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Symmetry Potential (E s ) [meV] 

FIG. 9: Eigenspectrum for a four donor cell as the symmetry 
potential (Es) is swept from -2 to 2 (meV) with an applied 
barrier potential (Eb) of 1.2meV 

putational state to the other. This is in exact anal- 
ogy with the technique of Rapid Adiabatic Passage for 
electromagnetically induced population transfer of atoms 
and molecules 25 . While this method has similar advan- 
tages to the adiabatic switching discussed by Lent and 
Tougaw^i, the method discussed here is entirely coher- 
ent and therefore doesn't rely on dissipation to ensure 
the ground state is always occupied. 

The effect of the barrier potential on the eigenspec- 
trum of a BDCA cell is shown in Fig. this time with 
an applied barrier potential of 1.2meV. The effect of this 
is to lower the barrier within each donor pair, delocalis- 
ing the electron and increasing the coupling between the 
donors. For a barrier potential of I.2meV there is now an 
energy gap between the ground and first excited states 
at the point of zero symmetry potential. At this point 
the eigenstates include contributions from all four basis 
states, not just the computational states. 

When Es = S max and Eb = 1.2meV the computa- 
tional ground state population has been reduced to 80% 
(compared to without barrier gate induced coupling be- 
tween the donors). For example, when Es — — 2meV 
the lowest energy eigenstate is 0.39|TT) + 0.89|TB) + 
0.U\BT) + 0.17\BB) whereas for E s = 2meV it is 
0.17|TT) + 0.U\TB)+0.89\BT) +0.39\BB) when written 
out in the position basis. 

The energy gap between the ground and first excited 
state allows adiabatic evolution to be used to shift the 
population from one computational state to the other, 
provided that the adiabatic criteria is satisfied. The adi- 
abatic criteria can be stated as2£ 



\(e\H\e) - (g\H\g)\ 2 
where \g) and |e) are the ground and first excited states 
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FIG. 10: The symmetry (Es) and splitting {Eg) potentials 
which are applied to the system to achieve adiabatic evolu- 
tion, where t p is the time over which the pulse is applied and 
the standard deviation of the pulses (a) is set to t p /6. B max 
and Smax are the maximum barrier and symmetry potentials 
and are used to control the amount of tunnelling and locali- 
sation respectively. 



FIG. 11: Eigenspectrum for a four donor cell as the pulse 
scheme given in Eq. 1121 and 1131 is applied to the barrier 
and symmetry gates respectively, where S max = 2meV and 
-Bmax = 1.2meV. 

Eq. (| X4p . to study the time dependence of the system 
including decoherence. The equation of motion is 



respectively. The pulse scheme given in Fig.^|can be ap- 
plied in order to achieve the energy level splitting while 
still ensuring the computational states are highly pop- 
ulated for readout. This involves applying a Gaussian 
pulse to the barrier gate while simultaneously switching 
the control gates from one polarity to the other. The fol- 
lowing functions were used for the barrier and symmetry 
gate potentials 



- n [H, P ] 



(14) 



E B (t) = S max exp 



Es{t) = S^erf 



2a 2 



(t - t p /2) 



V2cr 



(12) 



(13) 



where t p is the total time over which the pulse sequence 
is applied and a is the standard deviation of the pulse, 
which was set to a — t p /6. B max and >!? m ax are the max- 
imum barrier and symmetry potentials respectively. 

The resulting eigenspectrum is shown for this pulse se- 
quence in Fig. 1111 as a fraction of the pulse time t p . The 
degeneracy of the first two states is lifted but the compu- 
tational states are still strongly populated (> 99.999%) 
before and after the application of the pulse scheme as 
Eb ~ at t — and t — t p . 



VIII. TIME DEPENDANT BEHAVIOUR AND 
THE EFFECT OF DEPHASING 

As we are only considering a relatively small basis 
of states, we solve the density matrix master equation, 



where the Liouvillian [C[p]) describes the decoherence of 
the system. Integrating Eq. (|14|) in the limit of no deco- 
herence (C[p] = 0) gives the pure state of the system as it 
evolves over time. Fig. II 21 shows the state population for 
the pulse sequence given in section IV 111 over a pulse time 
t p = lOOps. The system is initially in \TB) and is then 
adiabatically switched to \BT) while only transiently oc- 
cupying the non-computational states. 

To determine the fidelity of transfer, we plot the fi- 
nal occupation probability of each state as a function of 
total pulse time (t p ) assuming we start in state \TB), 
Fig. ED Three distinct regions can be identified. For 
pulse times of less than O.lps, the pulse sequence is ap- 
plied too quickly for the system to evolve, which is to 
be expected as the pulse time is much less than the tun- 
nelling rates of the system. Pulse times of greater than 
20ps satisfy the adiabatic criteria, Eq. Ijlip. and the sys- 
tem moves smoothly from one computational state to the 
other with a fidelity of > 99.95%. Between these regions 
we see that after switching the non-computational states 
are occupied with varying probabilities. 

To study the effects of decoherence, we introduce 
£[p] ^ in Eq. I|14l) . As we have shown that the re- 
laxation rate (1/Ti) due to phonons is expected to be of 
the order of microseconds, we will only consider a phe- 
nomenological T 2 = 1/T 2 (pure dephasing) rate. We 
model this as a decay of the off-diagonal terms of the 
density matrix, 



£[p]=r 2 [p-diag(p)]. 



(15) 



8 



0.8 - 



°- 0.4 - 



0.2 - 



0.0 




20 40 60 80 100 

Time [ps] 

FIG. 12: (Color online) Population of states as a function 
of time showing complete population transfer from the state 
\TB) to state \BT) while only transiently populating the non- 
computational states (\BB) and \TT)). The time over which 
the pulse sequence is applied is t p — lOOps and the effects of 
decoherence are ignored. 
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FIG. 13: (Color online) Final population of states as a func- 
tion of total pulse time (t p ), ignoring the effects of decoher- 
ence. High fidelity transfer (> 99.95%) between computa- 
tional states is observed for pulse times greater than 20ps. 
For pulse times of less than O.lps the system does not have 
time to evolve from its initial state. Between these times, the 
non-computational states are partially occupied. 

Fig. El shows the probability of successful transfer from 
one computational state to the other as a function of 
total pulse time and dephasing time T^. The region of 
> 99% transfer is enclosed by the dotted contour line 
on the plot. The region in the top right corner corre- 
sponds to the system dephasing faster than it is being 
switched, resulting in a loss of fidelity. The final state 
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FIG. 14: (Color online) Probability of successful transfer as 
a function of both total pulse time and dephasing time. The 
region of > 99% successful transfer is enclosed by the dotted 
contour line in the bottom right corner. 

in this region is a uniform mixture of the computational 
and non-computational states. 

From Fig. [21 we see that even with finite dephasing, 
there is a window within which coherent transfer can still 
occur. For dephasing of 500ps and a total pulse time of 
20ps, a transfer probability of > 99% can be achieved. 



IX. SCALABILITY OF THE BURIED DONOR 
SCHEME 

The buried donor scheme can be scaled by adding more 
pairs of donors in a similar fashion to that used for the 
quantum-dot system to form a line of cells. The incoher- 
ent switching time for a line of cells is predicted to scale 
approximately linearly based on simulations of quantum- 
dot systems^. To compare this to the scaling of the co- 
herent scheme discussed in section lYlII we estimate the 
time to adiabatically switch a chain of buried-donor cells. 
As shown in Fig. this configuration involves a 'strip' 
barrier gate running the length of the chain and a pair of 
symmetry gates at one end. The switching of the chain is 
achieved by applying the same pulse sequence given ear- 
lier, Eq. H12fl and l|13fl . to coherently follow the adiabatic 
path from one computational state to the other. This is 
incorporated into an effective Hamiltonian of the form 

N N-l 

H cS = E s {t)a{ + E B {t) of + Jzz £ afaf +1 , (16) 

i=l i=l 

where N is the number of donor pairs. The interaction 
term Jzz is approximated using Eq. © where Jzz 

/ — E gs )/2 for a given cell size R. This model assumes 
that the entire chain forms a pure state throughout the 
transfer, which must be performed within the dephasing 
time of the system. In this case the minimum evolution 
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Fit Equation: 
B max =A*exp(-N/n)+BO 

BO 1.32 ±0.01 [meV] 
A 0.49 ±0.01 [meV] 
n 4.2 ±0.2 
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FIG. 15: Scaling behaviour for the maximum barrier cou- 
pling (.Bmax) which can be applied while still maintaining the 
minimum energy gap at the centre of the eigenspectrum (the 
degeneracy point). B ma x is calculated for increasing numbers 
of donors pairs and then fitted to an exponential function. 
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FIG. 16: Scaling behaviour of t a diab as function of number of 
donor pairs in a BDCA chain for various cell sizes (R). This 
gives an estimate for the scaling of the minimum allowable 
evolution time for high fidelity transfer. 



time will increase with the number of donor pairs and 
is controlled by the scaling behaviour of the energy gap 
between the ground state and the first excited state. The 
energy gap is limited by the height of the potential barrier 
the electron sees which is controlled by Eb- 

To simplify the analysis we will ignore decoherence and 
only consider situations where the minimum energy gap 
is positioned at the centre of the eigenspectrum (Fig. lTl"|) 
where Eg = 0. This puts a limit on the maximum bar- 
rier coupling (-B max ) which can be applied in order to 
introduce an energy gap. Fig. 1151 shows the maximum 
coupling (-Bmax) which still maintains the minimum en- 
ergy gap at the centre of the eigenspectrum. The scaling 
behaviour of the maximum allowable barrier coupling is 
fitted to an exponential decay, 

B max » 0.49 exp(-AT/4.2) + 1.32, (17) 

where N is the number of donor pairs and i? max is the 
maximum barrier coupling in meV. 

The exact behaviour of the system for large numbers 
of donors is computationally expensive to calculate. We 
can obtain an estimate for the minimum pulse time (t p ) 
which still provides high fidelity (> 99.95%) transfer by 
observing the scaling of the adiabatic time, 

6V25 max (ejof |ff) 

adiab ^ i£ gaP r ( ] 

with increasing number of donor pairs, where k is a scal- 
ing constant used to compare the minimum evolution 
time with the previous calculations. This equation is de- 
rived by assuming that the adiabatic criteria, Eq. (| 1 1 p . 
must be less than 1/k to achieve adiabatic evolution and 



noting that the time derivative of Eq. (|SJ| simplifies con- 
siderably at the degeneracy point (the centre of the eigen- 
spectrum, Fig. Ill|) . At this point the majority of the 
Hamiltonian is constant in time and the time derivative 
of the symmetry bias, Eq. (jl3[) . gives the numerical pref- 
actors of Eq. (jl8|l . We use k = 20ps as this is the required 
pulse time to achieve > 99.95% fidelity when switching a 
single cell comprised of 2 donor pairs with a cell size of 
R = 15nm. Fig. Ilblshows i a diab for up to 12 donor pairs 
(6 QDCA cells) calculated by diagonalising the effective 
Hamiltonian of the system, Eq. I|16|l . for several differ- 
ent cell sizes. As the cell size is increased, the coupling 
between pairs of donors is reduced and so the minimum 
switching time increases. 

While i a diab may underestimate the minimum allow- 
able evolution time, we expect the scaling behaviour of 
the system to be similar. The scaling behaviour is ap- 
proximately linear for these system parameters, though 
it will ultimately be restricted by the decoherence time 
of the system. This could be improved with a stronger 
bias field or by bringing the donors closer together. 

At present, more detailed simulations of logic struc- 
tures such as those shown in Fig.|21are required to demon- 
strate the viability of classical processing using BDCA. 
Based on simulations of incoherent adiabatic switching 26 , 
we expect the coherent control techniques discussed in 
section IVTT1 to be applicable to more complex logic struc- 
tures. 
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X. APPLICATION TO SOLID-STATE 
QUANTUM COMPUTING 

The use of BDCA have a number of applications for 
solid-state quantum computing, specifically silicon based 
architectures. Lines of BDCA cells could be used to ini- 
tialise charge based qubits and provide an interface be- 
tween the nano-scale features of the qubits themselves 
and the micro-scale control electronics. This would help 
to shield the qubits from the decoherence effects of the 
surrounding support electronics. The use of BDCA lines 
could also provide an on-chip pathway to transfer classi- 
cal information between nodes in a solid-state implemen- 
tation of a semi-classical or type-II quantum computer 28 . 

The adiabatic evolution discussed here could be ap- 
plied to many types of coherent systems based on 
quantum-dots, not just those based on buried donors. 
The adiabatic pathway allows for fast switching and high 
fidelity with minimal requirements on gate timing and ac- 
curacy and in this way is similar to Coherent Transfer by 
Adiabatic Passage 2 ^. The construction of an array of co- 
herent BDCA cells would also provide a demonstration 
of the scalability of semiconductor quantum computing 
schemes. 

The advantages of using adiabatic evolution (rather 
than coherent oscillations) to measure the decoherence 
properties of a charge qubit has been investigated by Bar- 
rett and MilburrjSS. The use of BDCA cells for this type 
of experiment has other advantages as well. A line of 
BDCA cells arranged with a series of sensitive electrom- 
eters along its length would enable the switching and de- 
coherence properties of the system to be more accurately 
measured. As there is more than one electron moving 
during a switching cycle, the use of correlated measure- 
ments between all of the electrometers (in the style of 
those used recently for detecting single electron transfer 
within a double-dot structur e)??'? 1 ) would provide a more 
accurate measurement than measuring the movement of 
just one electron. 



XI. CONCLUSION 

The use of buried dopants has been investigated as 
a possible implementation of quantum-dot cellular au- 
tomata, specifically for the case of phosphorous donors 
in silicon. The time and energy scales for this system 
were investigated and a model developed to describe the 
system evolution. 

For operating temperatures of less than 1A, the 
ground state occupation was found to be approximately 
100%. The incoherent switching time for a BDCA cell 
was found to be slow (of order microseconds) compared 
to other QDCA architectures due to the lack of a defined 
electron tunnelling pathway and poor coupling between 
cells. The use of coherent evolution to provide fast and 
controllable switching of BDCA cells was investigated for 
the case where quantum coherence can be maintained 



throughout the switching process. This was found to 
provide a fast and effective switching mechanism with a 
cell of side length 15nm having a coherent switching time 
of 20ps and a fidelity of greater than 99.95%. 

The effects of dephasing on this process was investi- 
gated and found to have minimal effect as long as the de- 
phasing time is approximately 10-100 times greater than 
the switching time. The scaling behaviour of the system 
was investigated for a simple line of cells and found to 
scale approximately linearly with the number of cells. 
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